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Too much noise kills entanglement. This is the main problem in its production and transmission. 
We use a handy approach to indicate noise resistance of entanglement of a bi-partite system described 
by d X d Hilbert space. Onr analysis uses a geometric approach based on the fact that if a scalar 
product of a vector s with a vector e is less than the square of the norm of e, then s yf e. We 
use such concepts for correlation tensors of separable and entangled states. As a general form 
correlation tensors for pairs of qudits, for d > 2, is very difficult to obtain, because one does not 
have a Bloch sphere for pure one qndit states, we use a simplified approach. The criterion reads: if 
the largest Schmidt eigenvalue of a correlation tensor is smaller than the square of its norm, then 
the state is entangled, this criterion is applied in the case of various types of noise admixtures to the 
initial (pure) state. These include white noise, colored noise, local depolarizing noise and amplitude 
damping noise. A broad set of numerical and analytical results is presented. As the other simple 
criterion for entanglement is violation of Bell’s inequalities, we also find critical noise parameters 
to violate specific family of Bell inequalities (CGLMP), for maximally entangled states. We give 
analytical forms of our results for d approaching infinity. 


I. INTRODUCTION 

Detection of entanglement is a fundamental problem in 
quantum information . Many theoretical works have been 
done to establish a criterion to identify multipartite en¬ 
tanglement Although not all of these works 

have direct impact on experiments, still they provide a 
mathematical background to understand the phenomena. 

An arbitrary density operator for two qudits can be 
written as 

d-1 

P = (Mo ®Mo+Y^ Ty M, 0 Mj (1) 

d^-i 

+ {d-l)Y^ {To^Mo 0Mi + T.oMi 0 Mq)) , 

i=l 

where Tij, Tgi, Tiq are components of the correlation ten¬ 
sor T, and (except f = 0) are the generalized Gell- 
Mann matrices (see Appendix [Al for details) with Mq = 

\J d{d-i) The components Ty are real and given by 

Ty = c(d)Tr[p(Mi 0 Mj)], where the coefficient c{d) de¬ 
pends on the dimensionality (see Appendix iBl) . 

In it has been shown that the correlation tensor 
could be utilized to construct entanglement indicators. 
The nonlinear character of the resulting entanglement 
indicator of make them more versatile than linear en¬ 
tanglement witnesses. The study of @ was mainly done 
for multi-qubit systems. Here, we use the geometric ap¬ 
proach of @ to detect entanglement for pairs of sys¬ 
tems, more complicated than qubits. In our approach we 
do not study the dynamics of noise, but we concentrate 
on properties of the state after evolving through a noisy 
channel. The white noise usually arises due to imper¬ 
fections in experimental set up, while a “colored” noise 
can appear when multi mrticle entanglement is produced 
via multiple emissions [ 3 ] . We also analyze the effects of 
different local noises, modeling random environment and 


dissipative process Q. 

We find that in some of the cases there exist significant 
differences in numerical values of the critical noise admix¬ 
ture parameters for having entanglement after a transfer 
via various noisy channels. A study of this kind was done 
for qubits in Ref. Q. Here, we study higher dimensions. 
Note that, a different approach to obtain the criterion of 
separab ility of entangled qutrits in noisy channel is given 
in Ref. [TJ. 

Another way to detect entanglement is to violate Bell 
inequalities [llj. In our paper, we discuss resistance to 
noise of Bell-type inequalities for arbitrary dimension d, 
e.g. Collins-Gisin-Linden-Massar-Popescu (CGLMP) in¬ 
equality M- 


II. SIMPLIFIED ENTANGLEMENT 
CONDITION 

Separable states can be described by a fully separable 
extended correlation tensor, T®®p = where 

yprod _ ^ r^k jg d^-dimensional gen¬ 

eralized Bloch vector of feth single qudit subsystem. To 
produce entanglement conditions we employ a geometri¬ 
cal approach based on correlation functions, which was 
introduced in Ref. The criterion has the following 
form: the state p described by its correlation tensor T is 
entangled, if 

max(f,f’P™'^) < (f,f’) = ||f||2. (2) 

(J^prod 

The scalar product may be defined in various ways, pro¬ 
vided it satisfies all required axioms. However the sim¬ 
plest choice is : 

d=-i 

{x,y)=J2x,,y,,. (3) 

i,i=l 
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In such a case, the maximum of the left-hand side 
of Eq. ([2]) can be obtained by the highest generalized 
Schmidt coefficient of the correlation tensor T as 

Emax = max (T,mi (g) m 2 ), (4) 

mi ,7722 


where fhk is a (d? — l)-dimensional generalized Bloch vec¬ 
tor which describes a quantum state of kth subsystem. 
Note that, for c? > 2 we do not have a Bloch sphere. Thus 
mfc’s represent the admissible Bloch vectors which repre¬ 
sent a physical state. One can also take the maximum of 
the left-hand side of (HI) as the square root of the highest 
eigenvalue of TT^, denoted by Tmax- However, note that 
Lmax > Tmaxj and this estimate gives a less robust crite¬ 
rion for entanglement. Still it is technically much easier, 
and thus we shall use it. 


III. BASIC AIMS 

We analyze entanglement of bipartite quantum states 
which are initially pure. For technical reasons we always 
put them in the Schmidt form: 


d-l 


IV') = 

2=0 


( 5 ) 


We study various noisy channels, i.e., white, product, col¬ 
ored, local depolarizing and amplitude damping. In some 
cases, one can define a critical parameter Ucrit for the 
noisy output, of the form p{v) = v\4’){4’\ -I- (1 — u)pnoise, 
such that for v > the state p{v) is entangled. 

The action of a channel on a density matrix of a two- 
qudit can also be written in terms of its operation ele¬ 
ments, more precisely p J2u,v ^uip) ® Ey{p)p{Eu{p) ® 
Ev{p))\ where Ei is an operation element of the channel 
parametrized by p, a certain strength parameter to be 
specified later. We define this parametrization in such a 
way, so that p = 1 defines no noise, just unitary(local) 
evolution, while p = 0 signifies the maximally noisy state 
leaving the channel. Hence, one can show presence of 
entanglement by satisfying the condition p > pE^. ent- 
Thus such a description is a generalization of the one 
involving Ucrit- We introduce a quantity 


iipiit,ENT) = min 


||f(p=l)|P ) 


( 6 ) 


to compare different noisy channels on the same footing. 
In order to check the entanglement condition (l2|), we put 
the correlation tensor of a state m in a Schmidt form 
itself. The general form of the tensor elements is given in 
Appendix [Cl Note that the Schmidt form of the state ([5]) 
does not imply a diagonal form of the corresponding ten¬ 
sor. Even for a two-qutrit state non-diagonal elements 
appear. This is not the case for qubits. Because of, the 


above reason, calculations of Tmax require the diagonal- 
ization of the tensor, which is a computationally difficult 
problem even for small d. 

For white, product, and colored noises we analytically 
compute the critical parameter for the maxi¬ 

mally entangled states (diagonal form of the correlation 
tensor)of a pair of qudits: 

d—1 

Kax) = 7^El**) (7) 

We present results for arbitrary d and for d —> 00 . 

We also analyze all non-maximally entangled pure ini¬ 
tial states of qutrits (d = 3), parametrized as follows: 

\'ijj^{a, (3)) = cosa|00) -I- sinasin/3|II) 

-I-sin a cos/3122). (8) 

The study is also done for lower Schmidt rank states, 
for d = 3, 4: 


d-l 

IV'rank-fc) = E (9) 

j=d-k 

where k < d and J2‘jZd-k = 1- 


IV. RESULTS FOR VARIOUS TYPES OF 
NOISE 


A. White noise 

White noise is represented by an admixture of the to¬ 
tally mixed state. A mixture of white noise and two-qudit 
pure state |'0) can be written as 

1 — V 

Pwhiteiv) + —^Id^^d- (10) 

Since white noise has no correlations, the correlation 
tensor T{v) for the state Pwhite('c) is given by T{v) = vT. 
In the case of the maximally entangled state, only diag¬ 
onal components of the tensor are nonzero and given by 
Tiiiy) = ±u/(d — 1) for i € {1, 2, ..., d^ — 1}. Therefore, 
we have Tmax = v/(d—l) and ||T(r;)|p = ?;^(d-|-l)/(d—1). 
Using the criterion ([2]), we can conclude that the state 
(HUD is entangled if u > (d -I- 1) This recovers the re¬ 
sult of Ref. [Q. In the limit of d —>■ oo, the critical value 
Kit, ENT tends to 0. 

For the case of two-qutrit non-maximally entangled 
states ([8|)we obtain Umax = v{d>+'/K)/Id and ||T(u)|p = 
v^{2 — A'/64), where AT = 31 -I- 12cos2a -|- 21cos4a -I- 
24 cos 4/3 sin^ a. In Fig. [1] we present the critical values 
^crit ENT- ^he pure state ([8]) is entangled for all a and /3 
except for a = 0, {a = = 0}, and {a = 7r/2,/3 = 

7r/2}. The lowest value of = 0.25 is obtained 

by the maximally entangled state |//’max)' 
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FIG. 1. The critical values Ucrit to detect entanglement 
for two-qutrit non-maximally entangled states ([S]) and white 
noise admixture. For v > the state is entangled. 

The lowest critical value, ent ~ 0.25, is obtained for the 
maximally entangled state IV'max)- 


FIG. 2. Gritical values Pcrit.BWT for entanglement detection 
are shown for two-qutrit non-maximally entangled state under 
local depolarizing noise. The lowest critical value, Pcrit.BiVT = 
0.5, is achieved for the maximally entangled state IV'max}- 


In Table El we present the numerical results for states 
dH). They are for two-qutrit states of the Schmidt form 
-^(|11)-|-122)), and of two-ququart states -^(|22)-h |33)) 

and + 122) -|- |33)). All this is compared with 

proper maximally entangled states. The analysis shows 
that the lower is the Schmidt rank of such states the more 
fragile is their entanglement. 


Therefore, we have Lmax = (l—r)^/(d—l) and ||T(r’)|p = 
(1 — r)^(d + l)/(d — 1). We define 1 — r = p, the critical 
value BATT local depolarized channel is obtained 
as 

Pait,ENT = (^ + 1) ^ ■ (14) 


1. White noise treated as local depolarizing noise 

A d-dimensional quantum state exiting a local depo¬ 
larizing channel is given by 

= ( 11 ) 

The use of such a formalism gives us a method of esti¬ 
mating how much noise (disturbance) is introduced per 
qudit. With probability 1 — r the initial state p still 
remain unaltered by the decoherence. To describe the 
local depolarizing effect on a two-qudit system, one can 
use the following transformation of the generalized Gell- 
Mann matrices: 


= ( 12 ) 

where /cS {l,2,...,(i^ — 1}, and make the replacement 
in Tij = c(d)Tr[p(Mi 0 Mj)]. 

For the maximally entangled state after such a 
transformation only diagonal components of the correla¬ 
tion tensor are nonvanishing, and they read 


T,.(r) =± 


(d-1) 


(13) 


We see that when d ^ oo, —>■ 0. Thus with 

increasing d the maximal entanglement is more and more 
resistant with respect to the depolarizing noise. 

Our study for arbitrary (initially) pure states of two 
qutrits (I5|) gives Tmax = p"^ {8 + v^)/16 and ||r(p)|P = 
p^(2 — Ar/64). The critical values Pent ent ploted 
in Fig. [21 It is entangled for all a and /3 except those 
related to product states: a = 0, {a = ^,/? = 0} and 
{a = ^,13= ^}. The lowest value of = 0.5 is 

obtained by the maximally entangled state ©. In Tab. m 
we give numerical results of the critical values p^nt ent 
for lower Schmidt rank states m for d = 3,4. We see 
that for a specific dimension, entangled states of lower 
Schmidt rank are less resistant to local depolarizing noise 
than entangled states of higher Schmidt rank. Thus the 
new analysis does not change the picture much. 

For local depolarizing channel ||'T(p)(j.j^ ent)\\'^ — 
and ||T(p = l)|p = for the maximally entangled 
state which yields '^(Fcrit ent) — This is identical 
to £; 7 VT detect entanglement for the maximally 
entangled state in a global depolarizing noisy channel 
(white noise). 


for i G {1, 2,..., — 1}. 
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FIG. 3. The critical values Ucrit.siVT to indicate entanglement 
for product noise. For a 0, {a —>■ 7r/2, /3 — ^ 0}, and {a —>■ 
it/2, P 7r/2}, the state (1151) is entangled for the critical 
visibility approaching 0.6039. 


B. Product noise 


where now the scalar product is defined by a positive 
semidefinite diagonal metric G 

d^-l 

(l,f)G= ^ (19) 

hi=i 

We consider a diagonal metric G with the following 
nonzero elements: Gu = 1 for i G {1,2,..., — 2} and 

Gd 2 _ij^ 2 _i = V. With such a metric, the left-hand side 
(fTSll reads L'z = u[l — v(d — 2)/(d — 1)1 and ||T'(u)||r 
ast;[l + t;(-6 + 6d-d2 + ,;(d_ 2 ) 2 )/(d-l) 2 ]. For d > 2, 
||T(z;)||q is always greater than implies that the 

state Pc('c) is always entangled except u = 0, i.e., even 
for an infinitesimally small v the state is entangled. 


D. Amplitude damping noise 

An amplitude damping channel can be described by a 
process of energy dissipation of a quantum system to en¬ 
vironment. The transition of excitation occurs between 
excited state and the ground state with a finite probabil¬ 
ity. No transition is allowed between excited states. Such 
an amplitude damping channel can be described by: 


We consider the following output states: 

Pprod(?^) = 'C|V’)('0I + (1 - v)pa 0 Pb, (15) 


d-1 

Pad(7’) Ek{r) 0 Em{r)p{Ek{r) 0 E^{r))\{20) 

k,m—0 


where pi = 1^® reduced density ma¬ 

trix of the ith subsystem. The critical values 
are presented in Fig. [3l The lowest critical value 
u^rit ent ~ obtained by the maximally entan¬ 

gled state IV’max) (for this state the product noise is re¬ 
duced to white noise) Note that = 1 for product 

states, whereas in the limits a —>■ 0, {a —>■ 7r/2,/3 —>• 0}, 
or {a —>■ tt/2,P — >■ 7r/2}, the state (I15p is entangled for 
the critical visibility approaching 0.6039. In other words, 
our geometrical method is very sensitive in the presence 
of product noise. 


C. Colored product noise 


In the special case of output states 

P^{V) = u|V'Lx)(^Lxl + (1 - 0 TTc, (16) 

with TTc = |c? — l)(d — 1|, the correlation tensor of the 
state (1161) has the following non-vanishing components: 


Tm('c) 


± 


d-l 


v{d—2) 

d-1 


for i G {1,..., — 2}, 

for i = d^ — 1. 


In order to get more efficient entanglement condition, we 
consider a generalized version of the criterion of Eq. @ 

max(f,fP™'')G< ||f||^, (18) 

'pprod 


where 


d-1 

Ao = |0)(0| + ^/^^^|^)(^|, 

2=1 


Ej = Vr \0){j\, j € {1,2, ...,d-1}. 


( 21 ) 


The transformation of the generalized Gell-Mann matri¬ 
ces is given by 

= tot M, £ , 

\ (1 — r)Mj otherwise. 


where j G (1, 2,..., d'^ — d}, fe G (2, 3,..., d}, the 
is the symmetric (or antisymmetric) Gell-Mann matrices 
(see Appendix lAl for details). Here, we do not consider 
the transformation rule for diagonal Gell-Mann matrices 
as they are redundant in later calculation. 

For the two-qudit maximally entangled state o , we 
have the following nonzero components of the correlation 
tensor 


Tu{r) 



( 1-0 


(d-1) 


2 


(d-1) 


for M, G 
otherwise, 


(23) 


where i G {l,2,...,d^ — d}. Although Tij{r) are non 
zero for i,j G {d^ — d, — d -I- 1, ...d^ — 1}, still they 
do not have any impact in our calculation to obtain 
the sufficient condition for entanglement. More pre¬ 
cisely, to calculate the critical value Pcrit for amplitude 
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damping noise with a specific combination of compo¬ 
nents of correlation tensor, we use the generalized cri¬ 
terion CHI) involving the diagonal metric; Gu = 1 for 
i G {1,2,..., — d} and other components are zero. As 

a result, we have Amax = (1 — f)/{d — 1) and ||T’(r)|p = 
(1—[2 -I- (d— 2)(1 — r)^] /(d—1). Thus, the entangle¬ 
ment criterion CHD reduces to (d — 2)p^ -|- 2p > 1, where 
p = 1 — r. We find that for d > 2 the state is entangled 
if the value p exceeds the critical value 

pint, ENT = (^2(d-2) ) 

where A = 2^/[3^(d — 2)] and B = [1 -|- (1 -I- 

The critical value goes to zero when d —>■ oo. 

It is worth noting that non-maximally entangled states 
for d = 3,4 in Eq. ® are more robust than the max¬ 
imally entangled state © against this type of noise. 
The corresponding non-maximally entangled states are 
l^(Tf’f)) for d = 3 and IV’nmax) = cos(0.853)|00) -b 
sin(0.853)(|ll) + \22) + |33)). The critical values 

Pcrit ENT ^-r® pfoted in Fig. ID In Table U we show nu¬ 
merical results of the critical values pf^,^ stvt for lower 
Schmidt rank states ([HI) for d = 3,4. It is found that the 
maximally entangled states o of lower Schmidt rank 
are less robust than the maximally entangled states of 
higher Schmidt rank. As a matter of fact the positive 
semidefinite diagonal metric G is not unique, one can ob¬ 
tain lower values for gjvT than the results presented 
in our manuscript with the optimal G. 

Also, we check critical parameter in ([HI) for the max¬ 
imally entangled state ©. Here = 

%^(2 + (d - 2){pi,,^p^py) and ||f (p = I)|p = 
^ which yields ^(Pcrit.BJVT) = Pcrit.EArT(2 + {d - 

‘^){piit,ENT?)^d-^. Hence, ^(Pcnt.^Arr) ^ 0, when d ^ 
oo. Due to our entanglement criterion we obtain lower 
value for ^(p^it ent) fo'^ ^ specific non-maximally entan¬ 
gled state Eq. ([5|) than the maximally entangled state d?]) 
while analyzing for amplitude damping noise. The cor¬ 
responding non-maximally entangled state is f)) 

for d = 3. For this specific state we ob¬ 
tain p^j.;^ = 0.4990. The critical values ^(p^it ent) 

are plotted in Fig. [51 A comparison between values 
of BATT ^{pirit ent) parameter, calculated for 
white or local depolarizing noise and amplitude damping 
noise, respectively is shown in Table [H]. 


V. TESTING ENTANGLEMENT WITH BELL 
INEQUALITIES 

In this section, we analyze noise resistance of violation 
of Bell-type inequalities for the bipartite quantum states 
©. To this end, we employ the CGLMP inequality in- 
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FIG. 4. We show the critical value PcTit,ENT of entanglement 
identification for non-maximally entangled state of qutrits for 
amplitude damping noise. The lowest value( 0.447(approxi- 
mate)) for pcrit, ent is obtained for the state |^(^, f ))(IH1) 
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FIG. 5. We show the critical value C(P?rit,EiVT) of en¬ 
tanglement identification for non-maximally entangled state 
of qutrits for amplitude damping noise. The lowest value 
(0.3560) for plrit,ENT is obtained for the state f ))® 


troduced in 


Ref. , which 


can be written as 



2k 

d- 1 


(Hfc-H_(fe+i)) <2, (25) 


where Bk = P{Ai = Bi + k) + P{Bi = A 2 + k A-1) + 
P{A2 = B 2 -l- fc) -l- = Ai k). 

We compare robustness of violation of the CGLMP in¬ 
equality, with respect to white noise, in case of the max¬ 
imally entangled state IV^^ax) ©i presented in [l^ with 
robustness of violation by the same state under different 
noisy channels. 
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noise 

d 

IV'max) 

IV’rank-2> 

IV’rank-s) 

l^nmax) 

d 

depolarizing noise 

amplitude damping noise 


2 

0.5773 

- 

- 

- 

2 

0.8410 

0.7071 

local depolarizing 

3 

0.5 

0.6546 

- 

- 

3 

0.8344 

0.7468 


4 

0.4472 

0.6794 

0.5283 

- 

4 

0.8310 

0.7647 


2 

0.5 

- 

- 

- 

5 

0.8290 

0.7750 

amplitude damping 

3 

0.4534 

0.8165 

- 

0.4465 

10 

0.8248 

0.7954 


4 

0.4239 

0.8660 

0.6124 

0.4074 

CX) 

0.8206 

0.8206 


TABLE 1. Summary of the critical values pfrit.BJVT to detect 
entanglement for the varions bipartite states under depolariz¬ 
ing and amplitude damping noises: the maximally entangled 
state I'i/'max)! the non-maximally entangled state |V'nmax)i a-nd 
the lower Schmidt rank states For amplitude damp¬ 

ing noise, the specific non-maximally entangled states present 
lower critical value Pcrit.BJVT than the maximally entangled 
ones. The specific states are presented below Eq. (HU. 


noise 

d 

IV'max) 

IV'i(ank-2) 

IV’rank-s) 

l^/'nmax) 


2 

0.3333 

- 

- 

- 

local depolarizing 

3 

0.25 

0.4285 

- 

- 


4 

0.2 

0.4615 

0.2790 



2 

0.5 

- 

- 

- 

amplitude damping 

3 

0.3888 

0.6667 

- 

0.3560 


4 

0.3256 

0.750 

0.3750 

- 


TABLE III. The critical parameter to indicate violation of 
CGLMP inequality, ptrit,LR are calculated for local depolariz¬ 
ing noise and amplitude damping noise for the bipartite max¬ 
imally entangled state in higher dimensional systems. 


that the lowest critical visibility for the CGLMP inequal¬ 
ity is achieved by non-maximally entangled states [l^ 
and studied in details (for d = 3) in Ref. [la. The 
lowest critical visibility r r;=0.6861')fRef. |l5j|') to 
violate CGLMP inequality in d = 3 is obtained for 
non-maximally entangled state |^^(1.0601, 0.7854)) (see 
([8])). For this specific non-maximally entangled state 
|'0^(1.0601,0.7854)), we obtain = 0.2883. 


1. White noise treated as local depolarizing channel 


TABLE 11. A summary of the critical values ^{pteit,ENT) for 
local depolarizing noises, and amplitude damping noise in de¬ 
tection of entanglement for ItAmax) and the lower Schmidt rank 
states IV'ink-fc) is given in a tabular form. As a mater of fact, 
^(Pcrit.BiVT) computed for each cases with local depolarizing 
noise is identical to obtained for the identical states 

with white noise. Also, for amplitude damping noise, the 
specific non-maximally entangled state presents lower critical 
values C(PCTit,BJVT) than the maximally entangled one. The 
specific state is Iz/'df, f)) for d = 3. 


A. White noise 

In this section we give an analytical solution for the 
critical visibility to violate the inequality [25] with the 
maximally entangled states and then consider the case 
for d —>• oo. 

Since white noise does not contribute to the violation 
of the CGLMP inequality (l25|l as stated in Ref. (l25l) . the 
quantum value for the state Pwhite('r') in Eq. (fTOl) is given 
by 


In this section we check the robustness of CGLMP in¬ 
equality in a local depolarizing channel with maximally 
entangled state IV'max)- 

Following the transformation rules for the Gell-Mann 
matrices shown in the state Iz/'max) after local depo¬ 
larizing channel can be written as 

Pdepoi(r) = (1 -IV^Lx)(V^Lxl + ^^^^2 (28) 

Since the noise part does not contribute to the Bell 
expression (l25|l with quantum mechanical observables, 
leads to the following condition for the critical parame¬ 
ters for the violation of (l25l) : 

Pir^t,LR = ( 29 ) 

where = (1 “ rt^t,LR)■ The results are given in 

Tab. mil However in terms of visibility the result 

is identical with (lT7l) . 


-^QM[Pwhite(l’)] — ■ydg^dz/;)]. (26) 

Therefore the critical parameters for the violation of the 
GGLMP inequality (l2^ are given by 

'^irit,LR = 2/7qm[|V')]- (27) 

and equal: « 0.7073, ^p 0.6962, 

<ru,LR ^ 0-6906 and « 0.6734 Note 


B. Colored noise 

In this section we check the violation of CGLMP in¬ 
equality by the maximally entangled state , \tpn-iax) CIIl 
mixed with colored noise 1161 Numerically we find that 
for d = 3 instead of choosing settings dm]) if we choose 
generalized SU{3) transformation matrix, then (1^51) is 
violated whenever ?; > 0. 
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C. Amplitude damping channel 


The same as in the previous cases, we examine the 
violation of the CGLMP inequality with state I'f/’max)- 
The state IV’max) after amplitude damping channel (|^ 
is given by 


PAD{r) = 

1 + r^(d — 1) 


/I \2 Cl—1 

|00)(00| + ^^ ^ \mm){kk\ 


m,k^0 


(l-r) 


r(l — r) 


d-l 


i-1 


|mm)(00| + \00){kk\ 

m^O k^O 


d-l 


d-l 


|TOO)(mO| + |0fc)(0fe| 

7717^0 k^O 


(30) 


We numerically check the critical parameters for d = 
2,3,4, 5,10. The results are given in Table Hill 

We see that the critical parameter increases with the 
dimension. It has a tendency to attain saturation for 
very large d. We give an analytical result for d tends 
to oo in Appendix ©• However, the lowest critical 
visibility (p^j.;j_^^=0.7316) to violate CGLMP inequal¬ 
ity in d = 3 is obtained for non-maximally entangled 
state |^/’^(0.8730,0.6449)) (see ([5])). For this specific non- 
maximally entangled state we obtain = 0.5, 

whereas ^(p^rit ent) — 0-4513. Also we check for 

I 


the specific non-maximally entangled state |'0^(ff i f)); 
which gives the lowest ent- We obtain = 

0.7429 for the state. We extend our calculation by check¬ 
ing p^j,;^ for the non-maximally state for 

which we obtain the lowest ent) considering all 

two-qutrit states. The value is obtained as 0.8640. 

Although for d —^ oo, p^^^ to violate (1^ in pres¬ 
ence of local depolarizing noise coincides (see Tab. [ml 
with the same parameter obtained in case of amplitude 
damping noise but in lower dimensions the variation in 
Pcrit LR icL case of amplitude damping noise differs with 
the one for local depolarizing noise. 


D. Fidelity of quantum states with critical noise 

In this section we give analytic expression of fidelity, 
i.e., the measure of closeness between two quantum 
states, e.g., the state after noisy channel with the pure 
state (here we only consider the maximally entangled 
state ([7])) of the system under investigation. 

Usually fidelity F for two quantum states is defined as 
a/ {(f>\p\4>), where ()> is a pure state whereas p is a mixed 
state. We choose © as (j), and p is parametrized by crit¬ 
ical parameters er ^r white noise and ^r for 
local depolarizing and amplitude damping noises. The 
general form of p’s are given in (fTOl) . (l28l) . and (l30l) for 
white , local depolarizing, and amplitude damping chan¬ 
nels, respectively. The corresponding critical fidelity is 
written as: 


jpd 

crit^LR 


crit,Li? 

y,d 


^CTit,LR \ i 


'W 


)^ 


/-/I „d A2 I ’'crit,l.B(2-'rcrit,I,il) a! 

lA — -I 

i((i+ id-mi - mm" +1)+ {d-mi - 


for white noise, 

for local depolarizing noise, (^f) 
for amplitude damping noise. 


We obtain identical value for lr fo^ white and local 
depolarizing channel but it differs with critical fidelity pa¬ 
rameter for amplitude damping noise in finite dimension. 
But, they coincide when d —> oo. To show comparison 
between critical values of fidelity parameter for different 
noises we write down our results for specific dimensions 
in a tabular form in Tab. Ml 


E. Werner gap 

In Tablel© we compare (j>irit,LR-Pirit,ENT) for local 
depolarizing and amplitude damping channels in case of 
bipartite maximally entangled states for d = 2,3,4,oo. 
As a matter of fact {ptrit.LR~ Ptrit,ENT) IIS'S same phys¬ 
ical interpretation as “Werner-gap”, which is defined as 
r’crit LR~''^crit ENT' ^ 1® apparent from the chart that, as 
we increase the dimension of the system, the difference 


d 

local depolarizing 

amplitude damping 

2 

0.8834 

0.8660 

3 

0.8544 

0.8397 

4 

0.8426 

0.8298 

OO 

0.8206 

0.8206 


TABLE IV. is calculated for white, local depolariz¬ 

ing and amplitude damping noise for the bipartite maximally 
entangled state. 


between pf^-^ er, and ^nt increases. We find that, 
when d ^ oo, p'^^^ ent f* for both noisy channels, 
whereas p'^i^R er increases with the increase of the di¬ 
mension of the system for amplitude damping noise and 
decreases with the increase of the dimension for depolar¬ 
izing noise. The overall effect reflects the fact that, if 






















d 

local depolarizing 

amplitude damping 

2 

0.2637 

0.2071 

3 

0.3344 

0.2934 

4 

0.3838 

0.3408 

CX) 

0.8206 

0.8206 


TABLE V. {pirit,LR — pirit.ENr) is Calculated for local depo¬ 
larizing noise and amplitude damping noise for the bipartite 
maximally entangled state. 


violation of CGLMP inequality is considered as an indi¬ 
cator to test the entanglement of the two-qudit maximally 
entangled state, then the maximally entangled state vio¬ 
lates the inequality for very high value of p'^rit lr com¬ 
parison to pirit,ENT- The discrepancy between Pcrit,LR^ 
and Pcrit ENT increases with the dimension. Therefore, 
the result indicates the known fact that not all states 
which are entangled violate specific Bell inequalities (in 
this case CGLMP ones). 

VI. FINAL REMARKS 

We analyze entanglement of bi-partite qudit systems 
using a special class of entanglement identifier defined 
with correlation tensors (l2|). To employ the criterion, 
in some cases we use operator-sum representations of 
the transformed states to describe the action of differ¬ 
ent noisy channels. We characterize the entanglement 
with a set of parameters, e.g., r'CTit.BArTi Pirit,ENT^ 
iiPcrit ent)- these parameters indicate identical be¬ 
havior of entanglement when d ^ oo. On the other hand 
we investigate violation of CGLMP inequality for the spe¬ 
cific states and also compare and ent- 

addition, we compute fidelity of states as a tool 

to quantify violation of local realism. At the end, it is 
worth to make few comments on the entanglement cri¬ 
terion with amplitude damping noise. May be due to 
the specific type of entanglement indicator, a specific 
non-maximally entangled state for d = 3 yields lower 
^{Pcrit ent)’ than the maximally entangled state of the 
respective dimension. It will also be interesting to in¬ 
vestigate if it can be possible to detect entanglement for 
lower values of ^(p^it ent)’ using either (12) or (fTOl) for 
different types of amplitude damping channels with dif¬ 
ferent types of Kraus operator representations other than 
the one used here. 
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Appendix A: Generalized Gell-Mann matrices 


The generalized Gell-Mann matrices in an arbitrary 
dimension d are the generators of the Lie algebra asso¬ 
ciated to the special unitary group SU(d). A recipe for 
the construction of the generalized Gell-Mann matrices 
is given in the next paragraph. 

Let Aj_fc denote a matrix with a 1 in the (j, k)- 
th entry and 0 elsewhere. This leads one to define 
three groups of matrices: (a) symmetric group; = 
Aj,fc + Afcj for 1 < j < /c < d. (b) antisymmet¬ 
ric group; = -i{Xj,k - Xkj) for 1 < j < fc < 

d. (c) diagonal group; “ 

ZA/+i_;+i) for 1 < Z < d — 1. The elements of three 
groups compose the generators of SU(d) group as Ad = 
{Mf, 2 , All® 3 ,..., Mi“ 2 , Aff^g,..., A/f,..., For an 

arbitrary dimension d, the total number of Gell-Mann 
matrices is d^ — 1. For simplicity, we re-index the ele¬ 
ments of the above set as Ad = 

The components satisfy the relations of tracelessness as 
Tr(AT,) = 0 and orthogonality as Tr{MiMj) = 25ij. For 
example, for d = 3 the set Ad is given by 


ALi = 


Afa = 


Ms = 


^0 1 o\ 
10 0 
»^o 0 oy 

^0 0 o\ 

0 0 1 , 
Vo 1 o) 

^0 0 -i\ 
0 0 0 
0 0 J 


, M2 = 


Mi = 


,Mfi = 


/o 

0 

A 

0 

0 

0 

VI 

0 

0; 

^0 

—i 

o\ 

i 

0 

0 


0 

0) 


0 

0 ^ 

0 

0 

—i 

Vo 

i 

OJ 



A 0 o\ 

0 

0 

M'j = 

II 

0 

1—1 

1 

0 

0 

0 


0 

0 

1 

. 0 

0 


(Al) 


The work is subsidized from funds for science for years 
2012-2015 approved for international co-financed project 
BRISQ2 by Polish Ministry of Science and Higher Edu¬ 
cation (MNiSW). AD was initially supported within the 


where the first three matrices, Mi, M 2 and M 3 , are sym¬ 
metric matrices, Mi,M^ and Mg are antisymmetric ma¬ 
trices, and M 7 and Mg are diagonal matrices. 
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Appendix B: Calculation of coefficient c{d) 

A single qudit density matrix is characterized by {(P — 
1 ) dimensional real vector as 

d^-i \ 

ld + k'Y^ PiMi j , (Bl) 

where pi = Tr(pMi) and we call p = [pi,..., pd' 2 _i) a 
Generalized Bloch vector, and fc is a normalization co¬ 
efficient that results in |^ = 1 for pure state. One can 
obtain the coefficient k = d{d — l )/2 as follows: for 
pure state, Tr(p^) = 1 /d + 2k'^\^^/d^ = 1. Therefore, 
the coefficient is given by fc = y/d{d — l)/ 2 . 

The element of correlation tensor for single qu¬ 
dit can be written as Tj = ^{d- i j [P^j\ — 

^(^'^i) Tr[/ 9 Mj]. Therefore the coefficient c{d) for bi¬ 
partite state is given by c{d) = 2 {d-i) • 



As previously mentioned in Sec. |TI1 the element of cor¬ 
relation tensor is given by T^- = c{d)Tr[p{Mi 0 Mj)], 
where Mi is the generalized Gell-Mann matrix presented 
in Appendix Here we do not consider the case where 
the elements of correlation tensor denote single particle 
correlation. We can divide the correlation tensor into two 
parts as follows; 


f = c{d) 


rpzonel Q 

Q ^zone2 


(Cl) 


where the is d{d — 1 ) x d{d — 1 ) matrix and its 

component is obtained by the Gell-Mann matrices in the 
symmetric and antisymmetric groups in Appendix [A) It 
has only diagonal elements as 


T^onei ^ ^2ajak for i G {1,..., d(d - 1)}, 


Appendix C: Correlation tensor for arbitrary 
bipartite state 

Here, we present a detailed form of elements of the cor¬ 
relation tensor for two-qudit states: |' 0 ) = X)n=o ctnlnn). 

I 


where j,k G { 0 , 1 ,..., d — 1 } and j < k. 

On one hand, the (d — 1) x (d — 1) matrix j^g^g 

not only off-diagonal components but also diagonal ones. 
They are obtained by the Gell-Mann matrices from the 
diagonal group. For the sake of simplicity, here, we re¬ 
index the elements of the matrix: 


/ 7 ~'zone 2 


7 ^zone 2 

J- fp-(^d-l),d'^-{d-2) 


'~pzone2 

-^d?-{d-l)4‘^-l 


\ 



T 


'Zone2 


a~'zone2 

-^d‘^-{d-2),d'^-{d-l) 


T 'zone2 

2,1 


^zone 2 \ 
-^ 1,^-1 


, 7 ^zone 2 
\ ^d'^-l4‘^-{d-l) 


7 ^zone 2 

^d'^- 14 '^-l 


/ 


7 ^zone2 
\ ^d-l,l 


7 ^zone 2 

^d-l,d-l 


/ 

(C2) 


Then, each element is given by 


outcomes with the following eigenvectors: 


T 


izone2 


i{i -I- 1 ) 




\n—0 


Tj 


izone2 


Vu(* + l)(j + l) V^o 

where i,j G {1,..., d — 1} and i < j. 




la. 


^zone2 _ ^zone2 

-^ji ~ ^ij ’ 


(C3) 


Appendix D: lQ^{d) for amplitude damping noise 

For amplitude damping noise, we choose the identical 
two measurement settings per p art y used to show vio¬ 
lation of CGLMP inequality [13. fig. The measurement 
operators Ag for Alice and Bt for Bob have d possible 


|^)a,s 

\b)B,t 


1 

Vd 


d-l 


j=0 


1 

7^ 


d-l 


3=0 


(Dl) 


where w = exp(27ri/d), a,b G {0,1,--- ,d — 1} and 
q;i = 0 ,Qf 2 = 1/2,/3i = 1/4, and P 2 = —1/4. This 
can be realized as follows: each one of the two parties 
first performs unitary operation on their subsystems with 

, , , r • 2m7r .2Tn(ci~l)7r ^ 

only diagonal entries - 3 where 

niG {as,Pt}- Then one of them performs a discrete 
Fourier transform Uft ( the matrix element is written 
as [UFT)pq = other one applies conjugate 

Uprp and at last they measure in their original bases. Ex¬ 
perimentally it can be realized by using phase shifter and 
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multi-port beam splitter on each side. Here, we denote 
Pq^{As = a, Bt = b) as the joint probability that Alice 
obtain the outcome a from the measurement As and Bob 
obtain the outcome b by measuring Bt for a given state 
Pad{t) in Eq. (l30l) . Then, the joint probability reads 


P^^iAs=a,Bt = b) 


1 — (d — l)(r — 2 )r ^ (1 — r) 


X 


d-1 

(l-r) uj^^-vh 

x,y^0 



(D2) 


where 'y = a — b +Ug + /3t- After summation, Eq. (ID2II is 
given by 


Pqm{^s — a, Bt — b) — 
(l-r) 


1 — (d — l)(r — 2 )r ^ (1 — r) 


d3 d3 

gin 2 sin 


Sin 




M7) 


- 1 


. (D3) 


To apply this joint probability into the CGLMP inequal¬ 
ity in Eq. (l25l) . we follow the definition introduced in 

Ref. [13: 


P{As 


d-1 

Bt + k) = YPi^s=n,Bt 

n—0 


n + k mod d). 


(D4) 


Also due to the symmetrical structure, the probability 
can be computed using simplihed definition: 


Pqm{^s = Bt + k) = ^ Pqj^{As = n,Bt = n + k) 
=0 

oAD 
QM 


n—0 

= dPA^{As = k,Bt = 0). (D5) 


This leads to the following relation: 


— Bi + k) — Pqm{Pi — A^ + k -\ 
= Pqm{M =B2 + k) 


1) 


= ^qm(52=Ai+A:). (D 6 ) 

Finally, for amplitude damping noise the quantum value 
of CGLMP inequality in Eq. (1^51) reads 


Ld/ 2 J -1 


/QM(rf,0=4 ^ 1 - 


fc =0 


2k 

d -1 


[qk{r) - g_(fe+i)(r)] , 
(D7) 

where gfc(r) = Pqm = Bi + k) and [pj implies the 
integer part of p. 

Now, we analyze the value of Eq. (ID7D for d —?► oo. The 
probability gfc(r) reads 


= lidAAhiAl + A;!) X 


( 1 -r)- 


d 2 

sm2 / ^G-i)(fc+i/4) \ 


d 2 


sm 


A(2d-l){k+l/4) \ 
d 


sm 


2 ( ■n-(fc+l/4) \ 


sm 


' 'K{k+l/A) ^ 


- 1 


In the limit of d —>■ oo, the only first sine function of (j'fc(r) 
survives and other terms reduce to zero. Finally, we have 


jAD 

^QM 


(oo,r) = 


2(1 — r )2 


E 


1 


1 


= 0-r)- 


k=0 

. 32 X Catalan 

’-- 

\2 


ik + ir ik + ir\ 


~ 2.9698(1 -r)G (D 8 ) 

where Catalan’s constant is Catalan ~ 0.9159. 
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